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ABSTRACT 

The time evolution of barred structures is examined under the influence of the external forces exerted 
by a spherical halo and by prolate halos. In particular, galaxy disks are placed in the plane including 
the major axis of prolate halos, whose configuration is often found in cosmological simulations. TV-body 
disks in fixed external halo fields are simulated, so that bars are formed via dynamical instability. In 
the subsequent evolution, the bars in prolate halos dissolve gradually with time, while the bar pattern 
in a spherical halo remains almost unchanged to the end of the simulation. The decay times of the bars 
suggest that they can be destroyed in a time smaller than a Hubble time. Our results indicate that this 
dissolution process could occur in real barred galaxies, if they are surrounded by massive dark prolate 
halos, and the configuration adopted here is not unusual from the viewpoint of galaxy formation. For a 
prolate halo model, an additional simulation that is restricted to two-dimensional in-plane motions has 
also ended up with similar bar dissolution. This means that the vertical motions of disk stars do not 
play an essential role in the bar dissolution demonstrated here. 

Subject headings: galaxies: evolution — galaxies: halos — galaxies: kinematics and dynamics — 
galaxies: structure — methods: n-body simulations 



1. INTRODUCTION 

Observations of disk galaxies have shown that thei r rota- 
tion curves are oft en flat out to large distances (e.g., Rubjg 
et al. 1982| , 1985; [Kent 1987| ). This finding suggests that 
disk galaxies, regardless of whether they are barred or non- 
barred, are surrounded by massive dark halo. However, it 
is difficult to distinguish between barred and non-barred 
galaxies from edge-on views which can easily provide ro- 
tation curve data. Recent development of observations 
enables us to measure rotation curves of barred galaxies, 
even though they are not viewed fully edge-on. As a re- 
sult, it has been confirmed that many, if not all, barred 
galaxies show fiat rotation curve, which does suggest the 
existence of massive dark halo in b arred galaxies as well 
QKormendy 1983| |Peters et al. 199l| ). 

The shape of dark matter halo, on the other hand, 
remains uncertai n observation ally, although some indica- 
tion is obtained ( [Sackctt 1999| ). Cosmological simulations 
demonstrate that dark matter halos are inclined to be 



more frequently prolate th an oblate ( Dubinski & Carlberg 



1991; Warren et al. 1992). In addition, prolate halos can 



be supported from the standpoint of the longevity of galac- 
tic warps; Ideta et al. (2000) have shown that a warped 
structure persists for a long time in a prolate halo while it 
disappears quickly in an oblate halo. 

Recently, El-Zant & Hafiler (1998) have revealed that 
bars are likely to be destroyed owing to chaotic diffusion 
of bar-supporting orbits, if they are embedded in triaxial 
halos including prolate configurations. Their results in- 
dicate that such bar destruction is conceivable from the 
situation where a barred galaxy is placed in a prolate halo 
with the major axis of the halo being in the disk plane. Ac- 



cording to Warren et al.'s (1992) cosmological simulations, 
this configuration of a disk in a prolate-like halo often oc- 
curs. Thus, real barred galaxies might be surrounded by 
prolate halos, and could suffer forces caused by the elon- 
gated potential of those halos. Unfortunately, however, 
El-Zant & Hafiler (1998) treated all the components of a 
bar, disk and halo as fixed potentials, and analyzed the 
orbits of test stars. Consequently, it is unclear whether 
bars are in reality destroyed in such prolate halos. 

In this Letter, we demonstrate from A^-body simulation 
that a barred structure is, in fact, completely destroyed 
in a time smaller than a Hubble time if it is embedded 
in a prolate halo whose major axis is placed in the disk 
plane. In §2, we describe the initial setup of our models 
and the numerical method. Results are presented in §3. A 
discussion and conclusions are given in §4. 

2. MODELS AND METHOD 

We set up a configuration such that the mid-plane of a 
galaxy disk is put in the plane including the major axis 
of a prolate halo. For comparison, a spherical halo is also 
included in our models. As a first step, the halo is handled 
as a fixed external field. We will discuss the possible effects 
of a live halo on the evolution of a bar pattern in the last 
section. The disk is evolved forward in time to form a bar 
via dynamical instability. Thereafter, we examine how the 
barred structure is influenced in the prolate halo. Since we 
pay attention to mainly the effects of an elongated halo on 
a bar, we do not include a bulge component which acts 
to reduce t he amplitudes of the bars arising fr om the bar 
instability (Hozumi, Fujiwara, & Nishida 1987). 

The prolate halo models are constructed from an ax- 



isymmetric modification of Hernquist's models (Hcrnquist 



2 



1990). The density profile of the halos is represented by 
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where My, is the halo mass, a and c are the scale lengths 
along the major and minor axes, respectively, and 



y 2 + z 2 



If the ratio of c to a is smaller than unity, the halo is 
prolate. We choose c/a — 0.6 and 0.75 as prolate ha- 
los, and c/a = 1 as a spherical halo. The models with 
c/a = 0.6, 0.75, and 1 are named Models P060, P075, and 
SI 00, respectively. 



We ado pt an exponential surface-density profile ( Frcc- 
man 1970) with a n isothermal d ensity distribution in the 



vertical direction ( Spitzer 1942 ) for the disk models that 
are given by 
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where Md is the disk mass, Rd is the scale length, and 
Zd is the scale height. The disks are truncated radially 
at 15 Rd, and vertically at 2 Zd- The velocity distribu- 
tion of disk particles is realized by employing Hcrnquist's 
(1993) approach which is based on the moments of the col- 
lisionless Boltzmann equation. In setting up the particle 
velocities, the disk plane is placed to include the major 
axis of the prolat e halo. The typical Toomre Q parame- 
ter ( Toomre 1964 ) is chosen to be of order unity so as to 
induce the bar instability easily, and the value, Q = 1, is 
set at the radius corresponding to the sun in the physical 
units of the Galaxy described below. 

TABLE 1 

Parameters for Jv-body Simulations 



Model 


M d 






M h 


a 


c 


P060 


1.0 


1.0 


0.2 


9.0 


10.0 


6.0 


P075 


1.0 


1.0 


0.2 


9.0 


10.0 


7.5 


S100 


1.0 


1.0 


0.2 


9.0 


10.0 


10.0 




Fig. 1. — Circular speed showing the contribution from the halo 
(dashed line) and disk (dotted line) to the total (solid line) for Model 
S100. 

We employ a system of units such that the gravitational 
constant G = 1, the disk mass Md = 1, and the scale 
length Rd = 1. If these units are scaled to physical val- 
ues appropriate for the Milky Way i.e., Rd — 3.5 kpc 
and Md = 5.6 x 10 10 M , unit time and velocity are 
1.31 x 10 7 yr and 262 kms -1 , respectively. The scale 
height Zd is set to be 0.2, or in our units, 700 pc, which cor- 
responds to that of old stars in the Milky Way. The disk is 



represented by 100,000 particles of equal mass. Since flat 
rotation curves of spiral galaxies suggest that the mean 
halo-to-disk mass ratio at t he Holmberg radius is 1.0 (e.g., 
Carignan & Freeman 1985), we determine the halo mass 
so that the disk and halo masses within 5 Rd are equal to 
each other in the spherical halo model. The parameters 
of each model are presented in Table 1, and the rotation 
curve of Model S100 is shown in Figure |l|. 

We use a hierarchical tree algorithm ( |Barnes fc Hut 
1986 ) with an opening angle criterion, 9 = 0.75, be- 
ing adopted. We expand forces and potentials up to 
quadrupole terms in the tree code . Forces are softened 



with a cubic spline ( Hcrnquist fe Katz 1989] ), and the 
spline softening length is 0.04 Rd, or in other words, 
0.2 Zd- The equations of motion are integrated with a 
fixed timestep, At — 0.04, using a leapfrog method. For 
these choices of parameters, the total energy was conserved 
to better than 0.13% in all simulations. 

3. RESULTS 

Figure || shows the contour plots of face-on density pro- 
files at intervals of 120 time units for Models P060, P075, 
and S100. In all the models, bars were formed through 
the bar instability, and they were almost fully developed 
by t ~ 120. We can see from Figure || that in the sub- 
sequent evolution, the bar pattern was kept nearly un- 
changed to the end of the simulation for the spherical halo 
model (Model S100), while the bars continued to dissolve 
gradually with time for the prolate halo models (Models 
P060 and P075). To evaluate the change in bar shape, we 
derived the axis ratio of the bars from the moment of in- 
ertia tensor for disk particles included in 2.0 Rd at which 
the bars end roughly. Then, we have found that the axis 
ratio changed from ~ 0.71 at t = 150 to > 0.95 at t = 600 
for Model P060, and that it did from ~ CK68 at t = 150 to 
> 0.93 at t = 600 for Model P075. 

To quantify the deformation of the bars formed in each 
halo model, we calculated their amplitudes as follows. The 
bars are very thin as compared to the extent of the disks, so 
that we ignore the thickness of the bars by projecting the 
particle distributions to the mid-plane of the disks. Then, 
we expand the density and potential of the projected dis- 
tribution in a set of basis functions as 



P ( R ) = X! A 
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where each density-potential pair, /_t nm and $ n rm consti- 
tutes Aoki & lye's (1978) basis set given by 
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Here, P nm are the Legendre functions (n > rn), and the 
radial transformation 
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Fig. 2. — Time evolution of the face-on density contours for Models P060 (top row), P075 (middle row), and S100 (bottom row). In all 
the models, bars are almost fully developed by t=120. The bar in the spherical halo model (Model S100) continues to exist to the end of the 
simulation, while the bars in the prolate halo models (Models P060 and P075) are weakened gradually with time and are completely destroyed 
at the end. 



is used. Positive values of m correspond to the number of 
arms in spiral patterns. In the expansions shown above, 
the amplitude of the (n, m)-mode is calculated from the 
absolute value of the expansion coefficient, |A nm |. In par- 
ticular, we pay attention to the fastest growing bar mode 
with (n,m) = (2,2). 




Fig. 3. — Time evolution of the bar amplitude of the fastest grow- 
ing mode for Models P060 (thin solid line), P075 (dashed line), and 
S100 (thick solid line). Also shown is the time evolution of the bar 
amplitude for a two-dimensional analogue of Model P060 (dotted 
line). 

In Figure we present the time evolution of the bar 
amplitude, 1 3.22 ] ; for each halo model. This figure de- 
scribes the behavior of the barred structures seen in Fig- 
ure H quantitatively. Again, we find that the bars have 
grown up to their full amplitudes by t ~ 120. In particu- 
lar, it should be emphasized that the bar amplitude decays 
nearly exponentially with time ~ exp(— t/r) for the pro- 
late halo models while it remains practically constant to 
the end of the simulation for the spherical halo model. To 
estimate the decay times r, least-squares fits were applied 
to the data from t = 200 to t = 400 for Model P060, and 
to those from t = 200 to 600 for Model P075. We ob- 
tained r ~ 111 and r ~ 375 for Models P060 and P075, 
respectively. If these values are represented by the physical 



units appropriate for the Milky Way, the decay times cor- 
respond to ~ 1.5 x 10 9 yr and ~ 4.9 x 10 9 yr, respectively. 
These time intervals are relatively smaller than a Hubble 
time that is almost equivalent to the supposed ages of disk 
galaxies, even though it takes about twice as long as the 
decay times for the bar amplitude to decrease by an order 
of magnitude. Thus, this kind of bar dissolution addressed 
here could occur in real barred galaxies, provided that our 
models reflect physical realism to a reasonable degree. 

4. DISCUSSION AND CONCLUSIONS 

We have found that a prolate halo can destroy a bar in 
a time smaller than a Hubble time, if the major axis of 
the halo lies in the disk plane. It has been known that 



a bar can be destructed by central massive objects (Nor 
man, Scllwood, & Hasan 1996; Hozumi & Hcrnquist 1999), 



and also by cl ose encounter or m inor merger of nearby 
dwarf galaxies ( Athanassoula 1996| ). Our finding adds an- 
other situation where bar dissolution can occur even with- 
out recourse to massive central objects or nearby compan- 
ion galaxies. Thus, in considering the evolution of barred 
galaxies, we should take into account the shape of sur- 
rounding halos and their configurations to the disks that 
support the bars. 

The time scale of bar dissolution will depend on the pa- 
rameters of halos such as the mass, core radius, and axis 
ratio. Then, to investigate whether real barred galaxies 
can experience bar dissolution, we need to check out how 
realistic our models are. First, for a halo model, we have 
used Hernquist's models that have a density profile propor- 
tional to r _1 near the center. Such a central cuspy density 
profile might be a natural end-product of hierarchical clus- 
tering based on a standard cold dark matter scenario (e.g., 
Navarro, Frcnk, fc White 1996| , 1997; [Fukushige k Makinc 
1997 ). From this point of view, our adopted halo density 
profiles would not be peculiar. Next, the halo mass and 
core radius were chosen to follow the observed rotation 
curves in spiral galaxies such that the circular velocities 



4 



are almost constant out to large radii. Thus, the kine- 
matic structures of the systems adopted here would not 
differ substantially from those of real galaxies. Last, Model 
P075 (c/a — 0.75) also leads to bar destruction within a 
Hubble time. Since cosmological simulations demonstrate 
that the angular momentum vector of galactic disks is of- 
ten aligned to the minor axis of dark halos, the axis ratio 
of halos in the disk plane corresponds to the major-to- 
intermediate one. According to Figure 7a of Warren et 
al. (1992), the cumulative fraction of the dark halos that 
have a major-to-intermediate axis ratio smaller than 0.75 
amounts to about 70%. In addition, the observationally 
determined geometrical form of dark halos is shown to be 
(b/a) p J> 0.8 as the equatorial axis ratio in density, and 
(c/a) p = 0. 5 ± 0.2 as the vertical-to-equatorial axis ratio 
in density (Sackett 1999). These values indicate that ha- 
los tend to be prolate with the major axis being in the 
disk plane, and that the axis ratio of c/a — 0.75 employed 
in our simulation is not far from the current lowest limit 
of 0.8. Taking into consideration those aspects mentioned 
above, bar dissolution might be in progress in the real uni- 
verse. 

As a possible mechanism of bar destruction, Raha et 
al. (1991) pointed out the bar buckling instability that is 
known as the firehose instability. This instability origi- 
nates in the anisotropic velocity dispersion between the 
motions parallel to and those vertical to the disk. If this 
is the case for the bar destruction presented here, a bar 
should survive in an infinitesimally thin disk, irrespective 
of the surrounding halo shape. Then, we run an additional 
simulation in which the motions of stars are restricted to 
the plan e of the disk with a self-consistent field (SCF) 



method dClutton-Brock 1972j ; pernquist & . Ostriker 1992| ) 
using Aoki & lye's basis set (eqs. || and Jt]|). The model 
realized with 100,000 particles of equal mass is a two- 
dimensional analogue of Model P060. In the SCF code, 
the maximum number of radial expansion coefficients is 
taken to be 16, and that of azimuthal ones is chosen to be 
2 with only even values being used (m = 0, and 2). The 
resulting behavior of the bar amplitude defined by IA22I is 
indicated by the dotted line in Figure |[ Clearly, the bar 
amplitude decreases nearly exponentially with time, and 
its time evolution is quite similar to that in the correspond- 



ing three-dimensional simulation (Model P060). This fact 
implies that the buckling instability should not be the es- 
sential cause of the bar destruction found in our three- 
dimensional simulations. 

Another likely mechanism of bar destruction is the ef- 
fect of chaotic orbits. Recently, El-Zant & Hafiler (1998) 
have found that the orbits that support a bar are highly 
chaotic when the barred galaxy is embedded in a triax- 
ial halo whose minor axis is vertical to the disk plane. 
Then, they suggest that the bar may be destroyed owing 
to chaotic diffusion, although they do not mention the de- 
struction time scale. Norman et al. (1996) have revealed 
that a bar is destroyed within a few dynamical times after 
the contraction of a rigid mass component, whose mass 
exceeds some critical value, is completed to form a central 
massive object. As a result, they are led to a picture that 
the abrupt destruction of a bar results from the chaotic 
behavior of bar-supporting orbits. However, our results 
show the gradual erosion of bars, and so, the behavior 
of the bar amplitude is closer to that found by Hozumi 
& Hernquist (1999) who examined secular change in bar 
pattern caused by central massive black holes added to flat 
disks. Of course, chaotic orbits might lead to the gradual 
destruction of bars. Therefore, we will analyze our models 
in phase space to study stochasticity in stellar orbits in a 
separate paper. 

Our simulations are highly idealized in that halos were 
treated as external fixed potentials. If halos are made mo- 
bile, the pattern speed of a bar will qui ckly slow down 



owing to dynamical friction of the halos (Weinberg 1985 
Debattista & Sellwood 1998). It is conceivable that the 



time scale of bar destruction could be affected by the pat- 
tern speed of a bar. Thus, a wider variety of simulations 
with a greater degree of realism will be required to figure 
out the detailed processes of bar dissolution. 
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